QCD factorisation and flavour symmetries illustrated in Bd^s KK decays 
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We present a new analysis of Bd,s KK modes within the SM, relating them in a controUed way 
through S'f/(3)-flavour symmetry and QCD-improved factorisation. We propose a set of sum rules 
for Bd,a K° K^ observables. We determine Bs KK branching ratios and CP-asymmetries as 
functions of Adir{Bd — > K^K'^), with a good agreement with current experimental measurements of 
CDF. Finally, we predict the amount of [/-spin breaking between Bd tt^tt" and Bs K^ K~ . 
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The current data in B-physics suggests that -B^-decays 
agree well with SM predictions, while B^-decays remain 
poorly known and might be affected by New Physics. 
Within the Standard Model, the CKM mechanism cor- 
relates the electroweak part of these transitions, but 
quantitative predictions are difficult due to hadronic ef- 
fects. The latter can be estimated relying on the approx- 
imate S'[/(3)-flavour symmetry of QCD : information on 
hadronic effects, extracted from data in one channel, can 
be exploited in other channels related by flavour sym- 
metry, leading to more accurate predictions within the 
Standard Model. 

In addition to isospin symmetry, an interesting theo- 
retical tool is provided by ?7-spin symmetry, which relates 
d- and s-quarks. Indeed, this symmetry holds for long- 
and short-distances and does not suffer from electroweak 
corrections, making it a valuable instrument to analyse 
processes with significant penguins and thus a potential 
sensitivity to New Physics. However, due to the signifi- 
cant difference — md, [/-spin breaking corrections of 
order 30 % may occur, depending on the processes. 

As a first application of [/-spin, relations were obtained 
between Bd — > tt+tt^ and Bs — > K^K^ . This led to cor- 
relations among the observables in the two decays such as 
branching ratios and CP asymmetries [ll, 0| and to a pre- 
diction for Br{Bs ^ K+R-) = (35t^g) • 10"^ [3. These 
results helped to investigate the potential of such decays 
to discover New Physics [3, Q . Unfortunately, the accu- 
racy of the method is limited not only by the persistent 
discrepancy between Babar and Belle on Bd tt+tt^ 
CP asymmetries, but also by poorly known [/-spin cor- 
rections. In these analyses, the ratio of tree contribu- 
tions Rc — |r^_|_/T^-|-| was taken from QCD sum rules as 
1.76±0.17, updated to 1.52t° }^ ^. In addition, the ratio 

of penguin-to-tree ratios f = \{Pk±/Tkp)/(P^±/T^±)\ 
was assumed equal to 1 [1] or 1 ± 0.2 [J, Q in agree- 
ment with rough estimates within QCD factorisation 
(QCDF) 3. Recent updates on [/-spin methods were 
given during this workshop [Sj]. 

QCDF may complement flavour symmetries by a more 
accurate study of short-distance effects. However, this 
expansion in as and 1 /mb cannot predict some significant 
l/mf,-suppressed long-distance effects, which have to be 



estimated through models. Recently, it was proposed to 
combine QCDF and [/-spin in the decays mediated by 
penguin operators Bd K°K° and Bs K°K° ^. 

The SM amplitude for a B decaying into two mesons 
can be split into tree and penguin contributions p^ : 



A = A{Bg ^ MM) = Ai«)r«f + A(«)p«f 



(1) 



with C denoting the charge of the decay products, and 
the products of CKM factors Ap'^ = ^pfc^pV Using 
QCDF [111, 113]' one can perform a l/mf,-expansion of 
the amplitude. The tree and penguin contributions in 
Bs K+K- and Bs K^K° in QCDF are : 
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- OL^EW + + 2/3" — -PsEW + 2(^4EW 



'^4EW + /^a + 2/^4 ~ 2P3EW + 2^4EW (3) 



' 2^'iEW ~ Paew (4) 

2^4EW + 1^3 + 2/3| — -P^EW ~ 04EW (5) 
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{0)fKGF/V^- The superscripts identify the 
channel and the bar denotes quantities for decays with a 
spectator s-quark. The tree and penguin contributions 
T'^° and P'^^ for Bd K°K° have the same structure 
as eqs. (jH) and ([51), with unbarred a's and /3's recalling 
the different nature of the spectator d-quark. 

At NLO in Us, a's are linear combinations of vertex 
corrections, hard-spectator terms and penguin contrac- 
tions, whereas /3's are sums of annihilation contributions. 
The weights of the various contributions are expressed in 
terms of as and Wilson coefficients [12]. a's and /3's con- 
tain the two most significant terms in the expan- 
sion: the LO terms, dominated by short distances, and 
the NLO terms in 1 /rrib that include the potentially large 
long-distance corrections. The latter, parameterised in 
QCDF through quantities denoted Xh (in power correc- 
tions to the hard-scattering part of ai) and Xa (in the 
annihilation parameters Pi), are singled out since they 
may upset the quick convergence of the l/mi, expansion. 
The other 1/mb-suppressed contributions, dominated by 
short distances, are under control and small, i.e, leading 
to a 0(5 — 10%) error. In 0, we showed that comparing 



FIG. 1: Diagrams contributing to Ba ~^ K"K° (left) and FIG. 2: Diagram contributing to ^ K+K from a tree 
Bs K°K^ (right) related through i7-spin transformations. operator, without counterpart in Bd ^ K'^K^. 



Bfi- and -B^-decays into the same final states helps to can- 
cel the potentially large long-distance l/mfc-suppressed 
effects {Xa,h), yielding improved SM predictions. 



I. SUM RULES 

Let us start with the difference = T'^^ — P'^'^ is free 
from the troublesome NLO infrared-divergence (modelled 
by Xa,h) that may be enhanced numerically by the chiral 
factor = 2m\/mi,/ms from twist-3 distribution am- 
plitudes. Hard-scattering {Xjj) and annihilation {Xa) 
terms occur in both penguin and tree contributions, but 
remarkably they cancel in the short-distance difference: 

Ad = [al - + - PI + 2Pl - 2/3^ (6) 
= Akk X asCpCi X [G{ml/ml) - G(0)]/(4^iV,) 

neglecting (small) electroweak contributions. The func- 
tion G ~ Gk 



^ Gk combines one- loop integrals from 
the penguin terms P4 and Pq defined in Sec 2.4 in ref. p^ . 
The same cancellation of long-distance 1 /wh-corrections 
happens for Ag = T"" — P"'^. Taking into account the 
uncertainties coming from the QCDF inputs we get 
Ad = (1.09±0.43)-10-^ + i(-3.02±0.97)-10-^GeV and 
As = (1.03 ± 0.41) • 10^^ + i(-2.85 ± 0.93) • 10-^GeV. 

These two theoretical quantities can be related to ob- 
servables, namely the corresponding branching ratio and 
coefficients of the time-dependent CP-asymmetry: 



T{Bd{t) ^ K^K°) - T{BS) A-°/?") 
T{Bd{t) ^ K^K^) + T{BS) ^ K^K^) 

_ AZ cos(AM • t) + <o„ sin(AM • t) 



(7) 



cosh(ArdV2) - Af sinh(Ardi/2) ' 

where we define [H: A^. = (|A|2 - |A|2)/(|A|2 + |A|2), 

+ =_ -(2e-^'^''AM)/(|A|2 -I- |A|2) and 0^ the 

phase of Bd — Bd mixing. A'^^ is unlikely to be measured 

due to the small width difference AF^, but it can be 

obtained from the other asymmetries by means of the 
relation K|2 + |A:^o^|2 + |A$^oj2^1, 

One can derive the following relation for Bd — > K^K^: 
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[X2 COS (pd + 2^3 sm (pd 



where Ld = Td^/U^ 
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4A'f|/(327rA'f|^) and: 



-1 = [|A('')p-t-|A('^)p-2|A('^)||AW|cos7]/n^ 

X2 = -[|A('^)|^ + |Ai'^)|^cos27-2|A('')||Ai'^)|cos7]/n^ 

X, = -[1-C0S7X lAl'^VlAi'^^ll/n, 



with n — 2 1 A, 
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A similar relation between 



As and Bs — > K^K'^ observables is obtained by replacing 
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and d s for all indices. 



These sum rules can be used either as a way to extract 
the SM value of one observable (say l^^irl) terms of the 
two others {BR"^ and Af^^^) and Ag, as a SM consistency 
test between BR^^ , \^dir \ ^'^'^ ^mix (and similarly for the 
B^ — > K^K^ observables), or as a way of determining 
CKM parameters [l^. These relations are free from the 
long-distance power-suppressed model-dependent quan- 
tities Xa and Xh that are a main error source in the 
direct computation of A'f^^ within QCDF. 



II. FLAVOUR SYMMETRIES AND QCDF 

Using U -spin symmetry, we can relate the two penguin- 
mediated decays Bd KoK n an d Bs KqKq, as exem- 
plified in fig. [1] (see also ref. [1^ in relation to ^ tttt) . 
[/-spin breaking should be much smaller here than usual: 
it does not affect final-state interaction since both decays 
involve the same outgoing state, and it shows up mainly 
in power-suppressed effects. This is confirmed by QCDF: 
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+6p^ + 26 PX - 6p^Ewl'^ - ^PXew) 



where we define the [/-spin breaking differences = 
— (id. for P). Apart from the factorisable ratio : 



/ = A^kk/AU = Mi^o"^^" (0)/[A/^P«'^-« (0)] 

which should be computed on the lattice, [/-spin breaking 
arises through 1/mb-suppressed contributions in which 
most long-distance contributions have cancelled out. 
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First, the hard-spectator scattering {Sa) probes the 
difference between B^- and i?s-distribution ampUtudes 
which is expected small, since the dynamics of the heavy- 
light meson in the limit nib — > oo should vary little from 
Bd and Bg. Second, the annihilation contributions (6(3) 
contain a [/-spin breaking part when the gluon is emit- 
ted from the light quark in the B^.s-meson (this effect 
from A\ and A| defined in [1^] is neglected in the QCDF 
model for annihilation terms). Taking the hadronic pa- 
rameters in [13, we obtain |P'*°/(/P'^") - 1| < 3% and 
|T^7(/T'^0) - 1| < 3%. 

Relations exist between Bd KqKq and Bg 
K'^K~ as well. A combination of J7-spin and isospin 
rotations leads from the penguin contribution in Bd ^ 
KoKq to that in B^ KqKq, then to Bs ^ K+ K' , 
up to electroweak corrections (it corresponds to fig. [1] 
up to replacing c? — > u in the right-hand diagram). On 
the other hand, there are no such relations between trees, 
since Bs — > K^K~ contains tree contributions (see fig. El) 
which have no counterpart in the penguin-mediated de- 
cay Bd KqKo- This is seen in QCDF as well: 
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Terms are ordered in decreasing size (in particular, curly 
brackets in T"^ should be tiny). From QCDF, we ob- 
tain the following bounds: |P"=^/(/P''°) - 1| < 2% and 
|T^±/(AJ^^ai) - 1 - r'^7(A^^ai)| < 4%. The lat- 
ter shows that flavour-symmetry breaking corrections are 
smaller than T'^° / [Aj^ j^ax) = 0(10%). Fortunately, 
is strongly CKM suppressed in Bg ~* K^K~ so that the 
uncertainty on its QCDF determination will affect the 
branching ratio and CP-asymmetries only marginally. 



III. SM PREDICTIONS FOR Bs KK DECAYS 

The dynamics of Bd K^K^ involves three hadronic 
real parameters (modulus of the tree, modulus of the 
penguin and relative phase) which we can pin down 
through three observables: BR'^° A'^^^ and A'^^^. Only 
BR'^^ = (0.96 ± 0.25) • 10"*^ ^ has been measured. 
However the direct asymmetry Afl^ should be observ- 
able fairly easily (for instance, A"^^^ = 0.19 ± 0.06 in 
QCDF) whereas the mixed asymmetry is likely small 
(^m°x = 0.05 ± 0.05 in QCDF). 

If only Af^^ becomes available, we have only 2 exper- 
imental constraints for 3 hadronic parameters. Then we 
may exploit a theoretically well-controlled QCDF con- 
straint to get T'^" and P'^" from BR'^°, Af^^ and the 



QCDF value of = T'^° — P'^^, free from infrared diver- 
gences and thus with little model dependence. This sys- 
tem yields two constraints in the complex plane {xp,yp) 
for P'^^ : a circular ring and a diagonal strip [9^], which 
can be satisfied only if \Afl^.\ < 0.2, and then yield two 
different solutions with opposite signs for Im P'"^, yield- 
ing two solutions for (P'"^, T'^°). 

From the measured value of the branching ratio for 
Bd — > K°K°, and choosing a particular value of the di- 
rect asymmetry A"^^^, we get the penguin and tree con- 
tributions as explained above. Then, the bounds in II 
yield the hadronic parameters in B^ — > KK decays up to 
small uncertainties. To be more conservative, we actually 
stretch the bounds in II relating Bd and Bg hadronic pa- 
rameters up to 5 % in order to account for well-behaved 
short-distance l/mfc-corrections not yet included. 

We obtain observables as functions of Af^l^, in Table [D 
In the case of the branching ratios, we have split the error 
in two parts. The first uncertainty comes from the QCDF 
estimates of Ad and ai, the theoretical constraints de- 
rived in II to relate Bd and Bg decays and the measure- 
ment of BR'^^ (this experimental uncertainty dominates 
the others). The second error stems from (factorisable) 
L/-spin breaking terms: / = 0.94 ± 0.2 (cf. 12]). 

Table [J corresponds only to the solution of the con- 
straints with Im P'^o > 0. But BR'^°, Af^^ and 
yield two different solutions for (r''°, P'^"), and thus 
for (r''±,P"±). Only one solution is physical, whereas 
the other stems from the non-linear nature of the con- 
straints. We can lift this ambiguity by exploiting a chan- 
nel related to Bg K^K~ through [/-spin, namely 
Bd TT+TT- [U, i, i, S, [i3- As explained in ref. [1], 
this allows us to reject one of the two solutions, and to 
compute for the other the J7-spin breaking parameters : 
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J- TTTT 



= 2.0±0.6 ^= =0.8±0.3 

^ 7^S± rTid± 

(11) 

The determination of BR"^^ is improved compared to 
the f7-spin extraction from Bd — *■ tt+tt" 0, 0, and in 
good agreement with the recent CDF measurement [l^: 



BR 



■s±\ 



I exp 



10" 



24.4 ± 1.4 ±4.6 



(12) 



IV. CONCLUSIONS 

We have combined experimental data, flavour symme- 
tries and QCDF to propose sum rules for Bd^g — > K'^K'^ 
observables and to give SM constraints on Bg KK 
in a controlled way. Tree (T"^") and penguin {P'^^) con- 
tributions to Bd K^K° can be determined by com- 
bining the currently available data with \T'^^ — P''"!, 
which can be accurately computed in QCDF because 
long-distance effects, seen as infrared divergences, can- 
cel in this difference, [/-spin suggests accurate relations 
between these hadronic parameters in Bd K^K^ and 
those in Bg K'^K^ . Actually, we expect similar long- 
distance effects since the K^K^ flnal state is invariant 
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Bi?"" X 10'^ 


Ail X 


X 102 


BR"^ X 10*^ 


X 102 


X 102 


Af> = -0.2 


18.4 ±6.5 ±3.6 


0.8 ±0.3 


-0.3 ±0.8 


21.9 ±7.9 ±4.3 


24.3 ± 18.4 


24.7 ±15.5 


= -0.1 


18.2 ±6.4 ±3.6 


0.4 ±0.3 


-0.7 ±0.7 


19.6 ±7.3 ±4.2 


35.7 ± 14.4 


7.7 ± 15.7 


Af =0 


18.1 ±6.3 ±3.6 


0±0.3 


-0.8 ±0.7 


17.8 ±6.0 ±3.7 


37.0 ± 12.3 


-9.3 ± 10.6 


= 0-1 


18.2 ±6.4 ±3.6 


-0.4 ±0.3 


-0.7 ±0.7 


16.4 ±5.7 ±3.3 


29.7 ± 19.9 


-26.3 ± 15.6 




18.4 ±6.5 ±3.6 


-0.8 ±0.3 


-0.3 ±0.8 


15.4 ±5.6 ±3.1 


6.8 ±28.9 


-40.2 ± 14.6 



TABLE I: Observables for Bs K'^K° and K^K as functions of the direct asymmetry Adir{Bd ~* K°K°) within the 

SM. We take Ai'*' = 0.0038 ■ e-'^, xi'^'' = -0.0094, Ai"' = 0.00088 ■ e"'^, A^"' = 0.04, and 7 = 62°, (/.^ = 47°, (Ps = -2° [H. 



under the d-s exchange. Short distances are also related 
since the two processes are mediated by penguin opera- 
tors through diagrams with the same topologies. J7-spin 
breaking arises only in a few places : factorisable correc- 
tions encoded in / = [M|^i^^=--^(0)]/[M|^F'»<'--^(0)], 
and non-factorisable corrections from weak annihilation 
and spectator scattering. 

Because of these expected tight relations, QCDF can 
be relied upon to assess t7-spin breaking between the 
two decays. Indeed, up to the factorisable factor /, pen- 
guin (as well as tree) contributions to both decays are 
numerically very close. Penguins in Bd — K'^K^ and 
Bs K~^K~ should have very close values as well, 
whereas no such relation exists for the (CKM-suppressed) 
tree contribution to the latter. 

These relations among hadronic parameters, inspired 
by L/-spin considerations and quantified within QCD 
factorisation, can be exploited to predict : Br{Bs — > 
K°K°) = (18±7±4±2)-10-6and Br{Bs K+K") = 
(20 ± 8 ± 4 ± 2) • 10^^, in very good agreement with 
the latest CDF measurement. The same method pro- 
vides significantly improved determinations of the U -spin 
breaking ratios ^ = 0.8 ± 0.3 and Rc = 2.0 ± 0.6. These 
results have been exploited to determine the impact of 
supersymmetric models on these decays [l9j . 

Our method merges ingredients from flavour symme- 
tries and QCD factorisation in order to improve the accu- 
racy of the predictions. Flavour symmetries must rely on 



global fudge factors typically of order 30%, without pro- 
viding hints where symmetry breaking is large or small 
(typically, is guesstimated, and not computed). QCD 
factorisation allows one to classify the contributions of 
the various operators of the effective Hamiltonian, but 
it suffers from a model-dependence in potentially large 
1/TOf, corrections. We use QCD factorisation to deter- 
mine the short-distance contributions, and we replace 
models for long distances by experimental pieces of infor- 
mation from other decays related by flavour symmetries. 
Both methods are exploited optimally to yield accurate 
predictions for Bs KK as functions of Bd observables. 

If sizeable NP effects occur, the SM correlations be- 
tween Bd and Bs decays exploited here should be bro- 
ken, leading to departure from our predictions. New re- 
sults on B ^ K form factors and on the Bd K^K° 
branching ratio and direct CP-asymmetry should lead to 
a significant improvement of the SM predictions in the 
Bs sector. The potential of other pairs of nonleptonic Bd 
and Bs decays remains to be investigated. 
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